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THE DYNAMICS OF OCEANIC FRONTS 
PART I THE GULF STREAM 

II 

Timothy W. Kao* 

Laboratory for Atmospheric Sciences (GLAS) 


ABSTRACT 

O 

The establishment and maintenance of the mean hydrographic properties of large scale density 
fronts in the upper ocean is considered. The dynamics is studied by posing an initial value problem 
starting with a near surface discharge of buoyant water with a prescribed density deficit into an, 
ambient stationary fluid of uniform density. The full time-dependent diffusion and Navier-Stbkes 
equations for a constant Coriolis p^ameter are used in this study. Scaling analysis reveals three 
independent length scales of the problem, namely a radius of deformation or inertial length scale, 

/j Ljj , a buoyancy length scale, h^, , and a diffus|ve length scale, h^. Two basic dimensionless para- 
meters are then formed from these length scales; the thermal (or more precisely, the densimetric) 
Rossby number, Ro= L^^/h^, and the Ekman number, E = The governing equations are 

then suitably scaled and the resulting normalized equations are shown to depend on E alone for^ 
problems of oceanic interest. Under this scaling, the solutions are similar for all Ro. It is also 
shown that 1/Ro is a measure of the frontal slope. The governing equations, in the form used in 
a previous paper by Kao, et al (1978b), are solved numerically and the scaling analysis is confirmed. 
The solution indicates that an equilibrium state is established. The front can then be rendered 
statioi|ary by a barotropic current from a larger scale along-front pressure gradient. In that quasi- ^ 
steady state, and for small valu^ps of E, the main thermocline and the inclined isopycnics fonning 
the front hav^ evolved, together with the along-front jet. Conservation of potential vorticity is 
also obtained in the light water pool. The surface jet exhibits anticyclonic shear in the light water 
pool and cyclonic shear across the front. The cross-front ageostrophic circulation is responsible 
for the maintenance of the front. It is also shown that horizontal diffusive effects are unimportant. 
Comparisons with known hydrographic features of the Gulf Stream are made, showing superb 
agreement. It is thus seen the mean Gulf Stream dynamics is indeed a solution of the Navier Stokes 
and diffusion equations. " 

For large values of E, it wUl be shown that another type of scaling is required. That result will 
be shown in a subsequent paper as Part II of this series, and is relevant to the study of density and 
current strocture on the east coast continental shelf of North America from Newfoundland to 
Chesapeake Bay, a region subject to forcing by freshwater river discharges. 
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THE DYNAMICS OF OCEANIC FRONTS 
PARTI THE GULF STREAM 


r. INTRODUCTION 

Large scale oceanic fronts are pseudo-pcnnanent features of the ocean, often representing the 
boundaries of water masses of different densities. The lighter water mass is bounded below by the 
base of the main oceanic thermocline. The isopycnics in the main thermocline are tilted upwards to 
form th^ inclined lateral boundary or the frbiit between the light water and the denser ambient 
water. There is therefore strong horizontal density change across the front. Obviously the horizontal 
scale of the frontal structure is much smaller than the larger scale of the water masses. Associated 
with the structure of the density anomaly is a geostrophic current along the front. The boundary 
currents on the western part of ocean basins are in fact such currents. A prominent example is the 
Gulf Stream in the North Atlantic. The lighter water in that case is provided by the Sargasso Sea. 
The inclined isopycnics of the density anomaly, i.e. the front, is sometimes referred to as the 
“wall” of the Gulf Stream. The isopycnics tilt upwards from the main thermocline at a slope 

•.V 

which averages approximately 1 :100. Clearly, it is the pressure gradient due to the buoyancy of 
the lighter water that is the immediate cause of the Gulf Stream. It was also observed (see Stommel, 

1 966) that the potential vorticity is a constant in the light water. Eased on these two concepts, 

' 1,7 

Stommel (see 1966) gave an inviscid, steady-state inertial theory of the Gulf Stream with a two- 
layer model. This simple model gave a good estimate for the upper ocean transport of the Gulf 
Stream but failed to produce the shape of tlie surface jet and over-estimated the maximum velocity 
of the current. To be sure, the model contained no information on the structure of the Gulf Stream. 
Later elaborations of the steady-state inertial theory by Chamey (1955) and Morgan (1956) 

incorporated the effects of the change of the Coriolis parameter, f, with latitude and the Sverdrup 

0 

transport from the interior ocean basin. 

In the present paper we investigate the dynamics of establishment and maintenance of the 
structure of oceanic fronts such as the Gulf Stream. The formation and maintenance of the Sargasso 
Sea is not addressed in that it is a part of the larger scale dynamics. The dynamics is studied by 
posing an initial value problem, starting with a near surface flow of light water arising from the 


1 


interior ocean basin of per unit length Into a stationary ambient denser nutd of unifonu density 
Pq. When equilibrium is establisliod the frontal structure advances steadily, except for a weak 
inertial period puliation. A barotropic current resulting from a larger scale along-front pressure 
gradient can be superposed to render the front stationary. This will concomittantly render the 
total vertically integrated cross-front transport to near 7.cro. It will be shown that the quasi-steady 
density and current structure of the front, obtained as solutions of the ftill Navie^$tokes and 
diffiision equations, possess a certain similarity property when proper scalings are used, The scaling 
analysis is at tlie heart of the present paper, Comparisons of the present results with the known 
hydrogtjaphic properties of the Gulf Stream are made. Excellent agreement is found. The results 
therefore elucidate the dynamics of the inean hydrography of the Gulf Stream. 

The inllow of buoyant surface water into au ambient stationary fluid has been subject to 
previous investigations by Kao. et .il ( l‘)77, 1 978a, b), hereafter referred to as I, II, and HI. We 
shall first summarize Uic physical understanding that we Ivave gained in these studies. 

(i) A surface density current whose front advances with a constant speed Ujj is established 
in time ~Ti (T| < I/O. The driving^^l'oree due to buoyancy is balanced at tliis stage by txictional 
dnig. The detailed distribution of the Isopycnics in the current depends only on the Reynolds 
number Re. (Re « wiiere v Is the kinematic viscosity). 

(ii) For times less than 1/f the effect of tlie eartli's rotation is largely neglegible. Thus in 
ostiun-inc plumes whose time of persistence is less than half of a semi-diurnal tidal period, the fron- 
tal structure is well-represeuted by the results without the earth’s rotation as given in I and in 
Garvinc (1974) for the ConnecEicut River plume reported by Garvine and Monk (1974). 

(iii) For longer times, the effect of the earth's rotation becomes increasingly more importiqjt. 
It becomes the dominant factor in time 2rr/f as sltovvn in III. The development of the along- 
front current calls into action tlie Coriolis Ibree which gradually relieves the frictional force from ^ 
its role in the balance of forces against the driving force of buoyancy. The forward speed of tlie 
front is greatly diminished and tlie stnictuix' tends toward a state ofgcostrophie balance. 

This is as far as we went In our previous st udies. Rigorous scaling considerations and detailed 
accounting of the relevance of the model to oceanographic events were not undertaken. In the [/ 
present paper, we first give the proper scaling laws. It Is found that tli4 problem is governed by a 


themial (or densimetric) Rossby number Ro and an Ekman number E. Furthermore, for sufficiently 
small values of E, a normalization can be found so that the solution is similar for all Ro and depends 

on B only, Results of the quasi-steady solution fora small B value are then given. 

A recent paper by Garvine (1979) on the steady-state frontal dynamics should be mentioned, 

Garvine used the momentum-integral technique, familiar in the engineering literature, by prescribing 
the density field and the shape of the velocity profiles. The prescribed cross-front velocity profile 
has a surface discharge towards the and a return flow at greater depth. The flow is assumed to 
be two-dimensional. Empirical interfacial fiiction and transport factors were assumed. The frontal 
shape and the along front velocity at the interface were then calculated. While the present study of 
frontal dynamics, within the framework of the fidl dynamical and thermodynamical equations, 
supports Garvine’s assumed surface discharge towards the front as a necessary ingradient for the 
maintenance of the front, the basis for Garvine’s assumption on the role of the intcrfacial stress is 
less apparent. The vertical Ekman number is a more natUTol pammetcr for measuring the effect of 
friction than Garvine’s scaling parameter Pr based on the interfacial stress coefficient. Indeed, in 
Garvine’s limit for small friction, Pr ■>«>, the conservation of potential vorticity in the light water 
side of the front was not obtained, whereas such a conservation is required by the dynainics, 

II. SCALING AND GENERAL CONSIDERATIONS 

We begin the scaling analysis by re-stating tlie initial value problem. A flow near the surface of 
Qg per unit length discharges into an ambient stationary fluid of density . The inflow has a den- 
sity deficit relative to the ambient of (Ap)^. The viscosity of the fluid is and the Schmidt num- 

ber is assumed to be unity. The local Coriolis parameter is denoted by f and g denotes the gravi- 
tational acceleration. The reduced gravity is g' with =^e,(Ap)Jp^, From and g' we form a 
buoyancy velocity scale , Uj = (g' and a buoyancy length scale h^ , h^ - (Q^/g')*'^ , The 
buoyancy time-scale (short time) is then Ti , Ti ~ h^,/U^j . From f we form an inertial time scale T, 

T - * /f, and an inertial length scale L^^ , L^, = U^/f. A diffusive length scale h^, is formed from 
and f, namely h^ = • Therefore there are three independent length scales to the problem 

Ljj, hg and h,,. These combine to form the two dimensionless parameters of the problem, i.e, the 
thermal or densimetric Rossby number Ro, 

« L 

Ro=r* 


and the Ekman number, £, 



(Note that Ko/E*Re, the Reynolds number, where Re* 


We let X* be the horkontal coH)rtlinale pointing opposite to the inflow direction and z* 
be the vertical co-ordinate measured upwards from the free surface. The flow is assumed to be in- 

dependent of the third co-ordinate y*. Let (u*, v*, w*) be the velocity components in the (x*,y*,z*) 

0 

directions. We now proceed to scale the problem. We let (t, » V (x'*'/Lg, z^/h^) t * t*f where t* 
is the dimensional time, and (ff, V, W)* (u*, v*, Ro w*)/0,, . The equation of continuity is then 


/ 
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dimensionless stream function is W, then 
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(1) 

( 2 ) 


We let the y* - component of the vorticity bo T*, i.c., t* = 3u*/3z* 3vv*/3x*. On normalizing f * 


by Uj/h^ we get, on using ( 2 ), 


i-«(Ro)-* 


ail .ait 

3^2 


(3) 


If we define the density anomaly as -y, 7 = (p - )/p^ and 7 ^ ® iAp\/p ^ , we can write tlie nomia- 

lized anomaly as, 


J-yhi, 


The governing equations are new: 


3y 3 a 


representing the diffusion equation, 


^ ^ 3^2 ^a^2 




3^“ 3p= 


(4) 

(5) 

( 6 ) 
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representing the Navier-Stokes equation governing the y component of the vorlictty, and[ 


( 7 ) 


representing the ^uation of motion for the V~ component of velocity. Equations (5), (6), (7) and 
(3) constitutes the complete set of equations, 

(The appropriate manner in which the non*rotating limit is to be viewed is outlined in Appendix B. 
It is seen that there the approach to steady»statc is not unifonn for Ro ». Thus for any non-zero 
rotation rate, f, however small, the steady-state solution for large times do«?is not approach the non- 
rotating limit. The non-rotating solution is however applicable for an intermediate time range 


Ti <t*<‘/0. 


/ 


For most problems of oceanographic ittterest Ro is large, say of the order of 10’ or larger. 
(We shall presently see that ‘ /Ro is in fact a representative slope of the frontal stnicture.) Thus we 
arrive at the result that the dynamics of tlie frontal structure is similar for all Ro and depends on 
the value of E only. This is so for all stages of the evolution of the frontal development. Indeed for 
Ro large the equations reduce to 


ay 3OT 

at at at? " an* 

a? 9(tfT) 

at a{ an an H an* 
at ai an an* 




an* 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 
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If we admit that the horizontal eddy coefficient greatei\than ( 1 ^^)^ (the vertical eddy co- 

efficient), it is clear that tlie similarity in Ro is still preserved unless (Vf)nKv^)y > (Ro)’ . This is 
highly unlikely to occur in tlie scale of motions that we are considering. 

We now return to show that * /Ro is in fact the representative slope of the front. This is easily 
shown to be a consequence of the thermal wind balance or Margules Law. Indeed from Margules 
Law for the limit of vanisliing viscosity, we have 






where s » frontal slope, and Av* is the across front gcostrophk velocity difference. Upon non- 
dimensionalizing the abp)%^ we have 

-'V / A7*Ros. 

We recall that the current was for short times t* < */f in force balance and was travelling at a steady 
speed U^ , This Is ttius an inertial current being turned by the earth’s rotation. At the new equilibrium 
state Av* ~ Uj or Av I . Thus 

S~*/Ro :7 

We Uierefore expect that the frontal slop will have a universal shape in the (t, v) co-ordinate in- 
dependent of Ro at any fixed t afterpuasl-gaostjrrtpliieequilibrium has been achieved, whehiquation 
(6) reduces to 

in the region where viscosity is unimportant, i.e. away from the surface Ekman layer. We antici- 
pate such a characterization of the slope to be valid for small values of E, E < 1 As E increases 
the region for which this characterization is valid decreases. At large values of £, £> 1 , the frontal 
structure is predominantly influenced by frictional effect. That case is of interest in the study of 
shelf-water density and circulation structures and will be investigated in Part II of the present 
series. It will be shown that another type of scaling law exists for that flow rdgime, 


It is also of interest to examine the vertical component of absolute vorticity, ai*, where w"* 


3v* 


= + f. Upon normalization by U^/h, we have 


9v 


On differenting Equation (7) with respect to ^ and 
the limit of vanishing viscosltyi 


! contfnwity equation, we get in 


For a region where ^ 0» |r easily seen to be also zero and the above equation reduces to 

9w.„9w. 3u 

+ M + o> 0 

3r 3^ 3| 

If such a region exists, and we will show in our initial value problem that it does, let us 
denote by D its dimensionless vertical extent. We then write the continuity equation as 

2!!+!JP.o 

d 3 3 3t D dt 

where ~ »^+ IT -r,.. Combining Equations ( 1 3) and ( 1 4) then yields the classical law of con- 
cii or of 

servation of potential vorticity 


( 12 ) 


(13) 


(14) 


.V' 
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s Q 


We thus see that the law of coniservaiion of potential vorticity is contained in the present formu* 
lotion. This law together with the thermal wind balance Is shown in the paper as a result of tlie 
initial value problem in the Un it of small viscosity. Thus the simple classical concept of the 
steady non-linear inertial the(|^y of the Gulf Stream as given in Stommel is embedded in the 
present framework, 


m. DISCUSSION OF RESULTS 


;v\ 


1 . Similarity and the Approach to Equilibrium 

Calculations were made using the formulation and numerical scheme given in III, to 
confirm the scaling laws, to investigate the approacli to quasi-steady equilibrium, and to give 
the density and current structure of the front. The relationship between the present normalized 
quantities and the corresponding quantities in the computational formulation is shown in 
Appendix A. In the actual numerical calculations, every term of the governing equations was, 
of coursCi included. 


Five cases were computed as summarized in Table I, The range of Ro is between 10 and 
210 and E is between 0.025 and 2.10. Cases 1 and 2 are identical in E but different in Ro and 

■H 


1 


0^ 




Table I 

Pirametrlc Regime of Computed Cases, 


Case No. 

Ro 

B 

1 (72C) 

26.06 

0.25 

2 (8IC) 

10,0 

0.25 

3 (80C) 

10.0 

0.025 

4 (82C) 

0 

26.06 

0.052 

S (78C) 

210.0 

2.10 


Cases 2 and 3 are Identical in ^o but different in E. Cases 4 and 5 differ in both Ro and E from 
every other case, 

The passage towards quasi-steadiness of the motion of the front is illustrated in Eigure i 
where the front speed Oj. is plotted against the dimensionless time r. The movement of the 
front undergoes a rapid deceleration in the first inertial day. Tlie adjustment towards geo- 
strophy is essentially completed at the end of this duration. This was already shown in III. 
Subsequent to this initial but primary adjustment, the speed of the front undergoes damped 
o.scillations of inertial frequency. The rate of damping depends on E, In Figure 1 the solid 
curve is for E « 0.025 while the dotted curve is for E * 0.25, i.e. a ten-fold increase in viscosity. 
Concomitant with tiic tioiizontal motion of the front, the frontal structure is deepening, as 
it must for mass conservation. The rate of deepening is also undergoing damped oscillations 
with inertial frequency but with opposite phase to the horizontal motion. Quasi-geostrophy 
is maintained in the frontal region away from the Ekman layer. At any time after the initial 
adjustment to geostrophy tlie shape of the inclined bounding isopycnic of the front, beneath 
the Ekman layer, is universal in the co-ordinates, defined by (5, t})/cc, where a is a con- 

stant of propotionality dependent only on the time, The bounding isopycnic of the front is 
defined here by 1/10 of the total density anomaly (Ap)Jp ^ . Tliis is shown in Figure 2 for the 
five cases listed in Table I at\ r ~ 7 with a = ir. It is seen that the shape is indeed universal. 


I. 
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Figure I. Front Speed 







For the larger values of E, tlu/Eknian layer is relailvely 4!hi0k and Is not exhibited In the Fiaure. 

(i 

For B ’<<] It the Ekinan layer tJiickness Is small. The dashed line Is the shape for 0.025. 

For that valuft of E, it is seen that the mean slope h/L of the bounding Isopycnic is equal to 
1.5 hg/Lj,, pr h/L » 1 .5 /Ro, where h and L are the vertical and horizontal pToJections of the 
inclined bounding isopycnic at any time. Thus h ** 1 .5 a h^ and L “aL^, Furthermore, this 
mean slope and the shape of the bounding isopycnic is preserved at a later time when quasi* 
steadiness is approached. The solid line in Figure 2 shows the shape at r ~ 30 (with 2sr). 

It is seen that the solid and dashed line arc almost coincident over its entirety. If wo denote the 
asymptotic values at equilibrium by an overbar, we then haveli« 1.5 ah,, and L*5L„ so that 
(H/L) remains approximately equal to 1.5/So for E 0.025. 


At any time after the initial adjustment to geostrophy, the structural features of the front, 
such as the distribution of the isopycnics forming tlic front, are found to depend only on E. 

For example, Figure 3 shows the surface manifestations of the density anomaly across the 
front for E « 0.025 and 0.25. They remain nearly identical for different values of Ro. Thus ' .c.' 
the scaling arguments advanced in the previous section are Indeed borne out by the results of 
the actual calculations, In addition, for small values of E, the features become not overly 
sensitive even to changes in E. For exaniple a change of E from 0.025 to 0.052 (case 4 of 
Table I) produces almost no change in the structure of the surface density anomaly. For larger 
values of E however, Figure 3 shows clearly that the compactness of the isopycnics at tlie front 
is governed by E. In general, increased vertical diffusion, as represented by the increase in E, 

causes a broadening of the isopycnics and a thickening of the Ekman layer, 

iJ 

The dynamics of establishment of equilibrium that emerges is really rather simple; namely 
after the initial primary adjustment to geostrophy, the shape of the frontal boundary in the 
(t«i »?«) ~ plane is universal away from the Ekman layer and is fonn-preserving in time. The 
front undergoes damped inertial period oscillations tending towards an equilibrium state when 
the deepening of the light water pool is restricted from below. The forward motion of the front 


CROSS FRONT SURFACE 
DENSITY ANOMALY 



Figure 3. Cross Front Surface Densiity Anomaly, 
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in the equilibrium state can ^fe brought to zero by adding a barotropic geostrojphic flow 

\ 

resulting from an along front pressure gradient of a larger scale motion. The structure of the 
front in the quasi-steady state depends only on E. Furthermore, for ^ 1, the structure is 

even insensitive to changes in the value of E. 

We now proceed to discuss the features nf ‘he quasi-steady solution for E ™ 0.025. In the 
next Section we will show that this case gives (iH ? pfincipal features of tlie mean hydrog- 
raphy of the Gulf Stream. 

2. The Quasi-Steady Solution and the Structural Features of the Front 

\\ 

The quasi-steady cross-sectional features of the frontal structure in the vertical plane is summa- 
rized and shown in Figure 4. The Ekman Number E is 0.025 and the solution is similar fcr all Ro*s. 
The outstanding feature is the compactness of the inclined isopycnics which form the front and the 
presence of a pool of light water of almost uniform density. The three dark solid lines are the 
isopycnics representing 10%, 50% and 90% of the total density anomaly between the incoming 
light water and the ambient fluid. (Hereafter we shall call them the 0.1, 0.5 and 0.9 isopycnics). 

To be more specific in tenninology we define the region bounded by the vertical at = 0 to the 

0. 1 isopycnic to be the frontal region. The part of the region of almost homogeneous light water 
bounded by the 0.9 isopycnic is the light water pool. The set of isopycnics between the 0. 1 to 0.9 
isopycnics, representing 80% of the total density anomaly, will be called the front. At “ 0, 
all the isopycnics become horizontal and the vertical extent bounded by the 0. 1 and 0,9 isopycnics 

V 

gives the tliickness of the main thermocline. We shall take the 0.9 isopycnic at ^g= 0 to be the 
base of the main thermocline. In Figure 4 the dashed curves are the isotachs of the along-front 
geostrophic current velocity?. The maximum difference in v across the front is 2.0 witli the maxi- 
mum occuring just inside the light water pool, The along front current forms a jet with cyclonic 
horizontal shear and strong vertical shear in the front, but anticyclonic horizontal shear and no 
vertical shear in the light water pool. Also shown in Figure 4, as represented by thin solid lines, 
are the instantaneous streamlines of the cross-front or lateral transport itt a frame of reference in 
which the front is stationary, The streamlines show surface convergence towards the front, down- 


13 


VERTICAL SECTION OF FRONTAL 

STRUCTURE 




- ISOPYCNICS 


V-ISOTACHS 

CROSS-FRONT STREAMLINES 


A 1*' 


••• D/Do (SEE TEXT) 

Figure 4. Vertical Section of Frontal Structure at Quasi-steady State, 





wolling and return flow at greater depths. It should be pointed out that In the quasi-steady state, 
the front Is still undergoing inertial perlorl oscillations, thougli of small amplitude. The relative 
streamline pattern in Figure 4 therefore changes over the period and so does the relative cross- 
front velocity (u- fff) even when the Isopycnics and along-front velocity structures are essentially 
unchanged. In any case, the relative cross-front velocities are typically 1/10 of Uj and the down* 
welling velocities are more than (Ko) * times smaller. Thus w* is 10"^ to 10’^ times the typical 
along-front surface geostrophic current. The features of the along-front jot and the density anomaly 
will now be examined in more detail. 

Figu re 5(a) shows the shape of the horizontal surface manifestation of the along-front jet, Tl»e 
jet is compact with a sharp peak at ~ 0.8 and is directed Northwards if the light water pool lies 
to the East of the front, It exhibits anticyclonic sheju* East of the front and cyclonic shear in the 

'.j 

front itself. The vertical profile of V'at location of the horizontal jet miaximum is shown in Figure 

5(b). The vertical shear in i^'is strongest under the jet maximum. We also observe from Figure 4 

9’v 

that inside the light pool v does not vary with deptli, so that ~is indeed zero there, We thcreloro 
expect the potential vorticity to bo conserved, which yields, together with geostrophy, a Stommel- 
type result that v decreases exponentially witli distance away from the front inside the light pool, 
and the depth of the light pool deepens expontentially. Indeed, we find 

and 

where is the dimensionless asymptotic depth of the light pool (D„ = above the main 
themiocline, i.e, « 0.82 according to Figure 4. The curve for D/D^ is shown by a dotted curve 
in Figure 4 and the curve for is shown by a short dashed curve in Figure 5(a). It should be 
noted that exponential dccajf ofV/vJ,j,jj and the exponential deepening of the light pool is generally 
faster than, tlvougli of tlie same order as, that of Stomniel (1966) for the same density anomaly aiid 

y/ 

the same asymptotic deptli. However, more than half of the total transport in Oie jet now takes 
place in the front itself, outside the potential vorticity conserving region of the light pool. The 
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Figure 5. Structure of Along-front Current: (a) Surface jet, (b) vertical profile at jet maximum, 






total transport is now found to be equal to 0,38 g'(D*)*/f compared with a value of 0,5(^ D*)’/f 
found by Stommel ( 1966). Although the transport is similar the velocity shear in the along front 
jet now changes from anticylonic shear in the light pool to cyclonic shear across the front. This 
feature, together with the conservation of potential vorticity inside the light pool, are benchmarks 
of success of the present formulation. 

Figure 6 shows two density anomaly profiles at station A-A and B>B of Figure 4. The density 
anomaly exhibits a deep homogeneous layer in the upper part followed by a well-defined thermo- 
cline. The proflle at A-A gives the thickness of the main thermocline whereas the profile B-B cuts 
across the sloping isopycnics of the front itself, 

Figure 7(a) and (b) show the surface expressions in density and sea-surface height anomalies 
respectively. The dimensionless sea-surface height anomaly, Ar;,, was found by integrating f 
vertically. If Ah, represents the dimensional rise of sea-surface height, then Ai?,- 
(Ah,/h) rp„/(Ap)g] , From the remote sensing, or satellite observational viewpoint, thesurface 
expressions, are rather important manifestations of the frontal structure, It is seen that the surface 
density anomaly varies more rapidly across the front than the sea-surface height anomaly. The 
latter varies over the entire width of the frontal region and contains information on the whole 
vertical cross-sectional structure of the front. The surface density anomaly on the other hand 
are easily influenced by events that take place in the surface Ekman layer, such as the effect of 
local wind and local heating and cooling. Thus altimetry data is generally a better indicator of 
the true location of the front than the infrared data on sea-surface temperature would indicate. 

The location of the steepest slope of the sea-surface height corresponds rather well with the location 
of the along-surface current jet maximum. Here again local surface wind effects in the surface 
Ekman layer has not been included. Local surface wind and heating effects are expected to alter 
the surface currents and surface temperature and salinity only in the surface Ekman layer. The 
importance of high quality altimeter data to examine the total depthwise frontal structure cannot 
be overemphasized. Indeed, given the maximum slope, the total rise height (Ah,)^,,^, and the 
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DENSITY ANOMALY PROFILES 
AT TWO LOCATIONS 


Figure 6. Density Anomaly Profiles at Two Cross-front Locations 
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Figure 7. Cross-front Surface Features; (a) density anomaly, (b) sea surface height anomaly 






crou-front width C of the i«a<Miiface height inomaly, the hydrography of the front it completely 
determined according to the present theory. 

fj^ 

IV. THE STRUCTURE OF THE GULF STREAM 

The results of § HI (2) are now applied to study the Gulf Stream. We confine ourself to the upper 
ocean of approximately 1,500 m or so. The top of the stratified deep ocean serves as a lower 
boundary of the upper ocean and imposes a limiting depth for the base of the main thermocline. 
Incorporation of this deep layer and the abyssmal current will be undertaken in a future paper. A 
few key references on the Gulf Stream should be mentioned: the book by Stommel (1966) and the 
references cited therein, the article by Fuglister (1963), the monograph by Worthington (1976), 
and a paper on numerjca! simulation by Semtnsr and Mints (1977). 

We may approach the application in the following way, We assume an imposed near surface 
transport towards the front, i.e. we assume a Q^, of given density deficit (Ap)Jp^ relative to the 
ambient water and examine the structure of the resultant density and flow fields in the quasi* 
steady state. From the structure, we obtain, in particular, the maximum surface jet velocity 

X\ (( 

the mean slope of the inclined isopycnics of the front, the width of the sea surface height anomaly 
L which Is the width of the frontal region, the maximum sea surface height anomaly (Ahj)^^,,^, and 
the Gulf Stream transport as functions of g', Q^, f and the depth of the upper ocean h. The for* 
mulas for these quantities as obtained from the results of § HI (2) are summarized in Table II. 

For example, we take as an initiai^^vhdition the inflow water to be Sargasso Sea water with a 
typical temperature of 20®C and 36S%o salinity and the ambient water to be the deep upper North 
Atlantic water of 5“C and 35.0%o salinity, so that Aa^ = 1.6 or (^p)Jp^= 1.6 X 10"^ and g' = 

1.5X1 0“* m/s^ . We take the magnitude of to be 50 m* /s and f to be 1 0“^ /s. The reference 
scales Uj, L^j, and h^ are then 0.91 m/s, 9.1 km, and 55 m respectively so that the densimetric 
Rossby number Ro = 165. To obtain the value of E, the Ekman number, we take the vertical eddy 
viscosity coefficient to be '^0(10) cm*/s. Then E ~ 0(10"’) to (KIO"’). Published estimates on 
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Table 11. 


/ List of formulas for Stream hydrography 


(( ' ‘ 

MAXIMUM GEOSTROPHIC SURFACE CURRENT, vj,, « 2 (g' 

(PI) 

FRONTAL SLOPE* 1.5 f (gV^/a Qj/a 

(F2) 

width of sea SURFACE HEIGHT ANOMALY, L r‘ (g')+2/a Q;H3 

1.5 

(F3) 

MAXIMUM SEA SURFACE HEIGHT ANOMALY « 0.7(gVg) K 

(F4) 

GULF STREAM TRANSPORT*P5|^«4-0.38£:HJ’ 

(F5) 


the magnitude of is subject to a range of variation, so that a precise specification is not justified, 

0 

nor is it necessary in the present framework since we have shown the results to be insensitive to 
changes in E for sufficiently iow values of E. On using the formulas listed in Table 11, with I *= 

1400 m, we get the maximum along-front Gulf Stream velocity v*,,^ to be 1.82 m/s, the mean slope 
of the inclined isopycnics to be 1 : 1 10, the width of the sea surface anomaly L to be 1 54 km, the 
maximum sea surface anomaly to be 1,57 m and the upper ocean Gulf Stream transport 

to be 33 X 10® m®/a. The details of the structural feature are shown in the similarity plot of 
Figure 4. From this we see that the width of the Gulf Stream itself, i.e. the width of the surface 
jet, is 108 km and the maximum depth of the 30 cm/s V- isotach is 750 m from the surface, 

The thickness of the main themiocline under the light water pool is 650 m. 

Table III shows the various Gulf Stream quantities for different values of while keeping 
g', f and h the same as in the example. These quantities were calculated using the formulas given 
in Table II, It is seen that a 1 5 times increase in increases the slope and maximum jet velocity 
by (15)*^^ or 2.47 times and decreases the width of the jet and L by the same amount, The Gulf 
Stream transport is of course independent of as seen from formula (F5) and depends only on 
h, g' and f. This must be so since the transport is induced gcostrophically from the buoyancy de< 
rived pressure gradient. For any , the Gulf Stream transport is sensitive to changes in E and g'. 
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Tible III. 

Table of Computed Gulf Stream quantitiei for different Q, 



INITIAL 

GULFS! 

rREAMQUANTITIBS 

Q, 

RBFERENC 

QUANTITI 

:e 

ES 

!lo 

MEAN 
SLOPE OF 

(hh 

WIDTH OF 

SURFACE JET 

widtoopsba- 

SURFACE HEIGHT 

(m*/s) 

PI 


K 

(m) 


FRONT 

(km) 

anomaly 

C(km) 

10 

0.53 

5.3 

19 

279 

111 86 

1.08 

182 

260 

25 

0.72 

7,2 

35 

205 

lil36 

1.44 

136 

194 

50 ‘ 

0.90 

9.0 

55 

165 

hllO 

1.80 

107 

154 

75 

1.03 

10.3 

73 

142 

li95 

2.06 

94 

134 

100 

1.14 

11.4 

88 

129 


2.28 

85 

121 

125 

L23 

12.3 

102 

120 

1;80 

2.46 

79 

113 

150 

1.30 

13.0 

115 

113 

1:75 

2.60 

74 

105 


Thus, for comparison with Stommcl’s (1966) result, taking the depth of the homogeneous light 
pool to be 800 m and g' to be 2 X 10"* m/s* as used by Stommel ( 1 966), we get a value for the 
transport of about 50 X 10* m*/s, The values given for the Gulf Stream quantities in Table III 
are all accessible to the Stream. This then suggests wide variations of Q^. Indeed variations of 
by several times is reasonable from considering time-dependent variations of wind-driven contri- 
bution alone. In addition, there is possibly a significant i^iintribution from the Gulf Stream return 
flow as described by Worthington (1976). 

In the present framework, we now see that the Gulf Stream dynamics is maintained by the 
near surface shoreward transport of light Sargasso Sea water. The Gulf Stream is a natural geo- 
strophic response to the pressure gradient associated with this buoyant transport. However, at any 
typical section, say at Cape Fear, the lateral or cross-front forcing induces a geostrophic Gulf 
Stream transport which is fed by water from further upstream, i.e. from the South. The com- 
position of the water in the Gulf Stream therefore reflects its tropical orgin. Between two neigh- 

I 

boring sections, say between Woods Hole and Hatteras, the increase in Gulf Stream transport, if 















































any, must be contributed by entraining water between the two sections from the Sargasso Sea 
which is the parent pool of light water. This concept appears to be compatible with the tight 
re^jirculation pattern of the Gulf Stream proposed by Worthington (1976). Finally, it should be 

'\V 

pointed out that the formation of the Sargasso Sea warm lens is perhaps partly wind-driven as 
shown by Semtner and Mintz ( 1 977), but its study belongs to the study of larger scale events and 
is beyond the scope of the study of frontal dynamics of the present paper. 
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, APPENDIX A 

In tlio fonnuiation given in III, the region of integration has a finite depth d. For the sake of com- 
putational convenience, d is used as the reference length, QJA as the reference velocity and 
d*/Q^ as the reference time in that scheme. The dimensionless parameters are Ro *U/(fd), Re “ 
Ud/nj,, F*U/(gd)*/2 and » 10^/* F/^^. The conversions of the computer results to the present 


normalization is achieved by the following fonnulas; 

t * (lORor* f/'^ X (AI) 

n *10F,-2/3 j. (A2) 

U » (10)"^ F^2/3 (^3) 

w »10 RoI\“2/3w (A4) 

r « (Ro)^i t (AS) 

Ro » lOORolPjl^ (A6) 

E « 100Ro/(ReF//3) (A?) 

hjd« (10~>)F^2/3 (AS) 


A"1 
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APPENDIX B 


If we normalize the governing equations by the buoyancy velocity, length and time scales, the 
equations corresponding to (3), (5), (6) and (V) arc 


3f 



(Bl) 


(B2) 

' n, Ro a, at, Rc ’ ^ 

(B3) 


(B4) 


3 ® . 32 


where = K */\ , tj - t*/Ti , Wj “ w^/U^ v* s ^ + 


9^1 


, 2 ' 


In the non-rotating limit Ro = ««, (B4) Is uncoupled from the other equations, and is identically 
zero for an initial value problem starting from rest. The steady-state is approached for ti °° , 

T = 0. For Ro # oo however, the steady-state is approached for tj -»■ «> and r large. The approach 
to steady-state is thus not uniform for Ro -*• Hence the limit tj -► <» and Ko -+ «» are not inter- 
changeable. 
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